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ABSTRACT 

We study the X-ray cluster gas density distribution in hydrostatic equilibrium us- 
ing the universal temperature profile obtained from recent simulations involving only 
gravitational processes. If this temperature profile is an indicator of the influence of 
gravitational processes alone on the intracluster medium, then the comparison of var- 
ious X-ray parameters expected from this profile and the observed data would point 
towards any additional physics that may be required. We compare the entropy at 
0.1i?2O0 and i?5oo, the scaled entropy profile, the gas fraction at 0.3i?2oo and the gas 
fraction profile with recent observations and discuss the implications of this tempera- 
ture profile in light of these data. We find that the entropy imparted to the gas from 
gravitational processes alone is larger than previously thought. The entropy at -R500 
for rich clusters is consistent with data, whereas the entropy at 0.1i?2oo is still less 
than the observed values. We also find that the gas fraction in the inner region of clus- 
ters, expected from gravitational processes alone, is smaller than previously thought 
but larger than the observed data. It does show a trend with the emission-weighted 
temperature ((T)) as shown by data. We therefore find that the role of any addi- 
tional non-gravitational process influencing the physical state of ICM would have to 
be revised in light of these findings. 
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1 INTRODUCTION 

The formation of structures in the Universe is believed to be 
hierarchical, as primordial density fluctuations, amplified by 
gravity, collapse and merge to form progressively larger sys- 
tems. This hierarchical development leads to the prediction 
of self-similar scalings between systems of different masses 
and at different epochs. 

X-ray observations of the intra-cluster medium (ICM) 
provide an ideal probe to test these self-similar scalings. Ob- 
servations however present a picture that is at variance with 
predictions from these self-similar scalings. The X-ray lumi- 
nosity of low temperature (poor) clusters fall below the self- 
similar expectations, and the surface brightness profiles of 
these clusters are seen to be shallower than those of richer 
clusters. It is also instructive to view this in terms of the 
entropy of the intra-cluster gas, which in self-similar scaling 
should increase in a very simple scaling with the mean tem- 
perature of virialized systems, whereas observations show 
that gas in low temperature clusters have larger entropy 
than expected. Another probe of non-gravitational processes 
is the gas fraction by mass in these clusters. Observations 



show that the inner regions of poor clusters have less gas 
(compared to the total mass within those radii) than those 
of rich clusters. 

These differences between theoretical expectations and 
observations have led to the emergence of a number of theo- 
retical ideas to increase the entropy of the ICM by some non- 
gravitational process, especially in low temperature clusters, 
involving heat input from supernovae (Valageas & Silk 1999; 
Kravtsov & Yepes 2000; Wu et al. 2000), quasar outflows 
(Bower 1997; Loewenstein 2000; Voit & Bryan 2001; Nath & 
Roychowdhury 2001), gas cooling (Knight & Ponman 1997; 
Bryan 2000; Voit & Bryan 2001; Pearce et al. 2000; Muan- 
wong et al. 2001; Wu & Xue 2002a; Dave et al. 2002) and 
accretion shocks (Tozzi & Norman 2001; Dos Santos & Dore 
2002; Babul et al. 2002), though it appears that there are 
problems with many of them (Ponman et al. 2003 and refer- 
ences therein; hereafter PSF03). 

A recent result from a very high resolution numeri- 
cal simulation (Loken et al. 2002) offers a fresh look at the 
scaling laws expected from gravitational interactions. Lo- 
ken et al. (2002) found a universal temperature proflle for 
the ICM, when scaled to the emission weighted tempera- 
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ture, from the assumptions of pure gravitational evolution 
of the cluster gas, with no input from any additional non- 
gravitational processes of heating or cooling. Interestingly, 
this temperature profile is a good match to the observed 
universal temperature profile, leaving aside the very central 
regions (r < O.lRyir) (De Grandi & Molendi 2002). 

The interesting aspect of this temperature profile ob- 
tained from this simulation is that it has a core within which 
the profile is fiatter than in the outer region. It has been ear- 
lier noted by Eke ct al. (1998) that interactions between dark 
matter and baryons imparts some entropy to the gas pro- 
viding it with a core (and making it deviate in the central 
regions from the self-similar expectations in which gas den- 
sity is proportional to the dark matter density). The entropy 
imparted to the gas was, however, thought to be important 
only in the very central regions (r <C 0.1ii„ir),and observa- 
tions of entropy at ~ O.lRvir was thought to refiect other, 
non-gravitational influences, if there wore any. 

It is then important to determine the gas density profile 
corresponding to this temperature profile, assuming hydro- 
static equilibrium, and determine the level of entropy ob- 
tained only from gravitational interactions, without the aid 
of non-gravitational processes. In this paper, we study var- 
ious implications of this temperature profile, especially in 
light of the entropy problem mentioned earlier. 

As we describe below, we find that the entropy imparted 
to the gas by gravitational interactions alone is larger than 
previously thought for poor clusters, even at radii exceeding 
0.17?2oo, and this is the main result of the present paper. We 
begin with the description of the background dark matter 
density profile that is assumed (§2). We then determine the 
gas density profile assuming hydrostatic equilibrium in §3 
and discuss various aspects of our result, including compar- 
ison with relevant data in §4. We end with a brief discussion 
of our results in §5. 

We assume throughout the paper Oa ~ 0.71, flo = 
0.29, Oi, = 0.047 and h = 0.71 which are the best fit pa- 
rameter from WMAP (Spergel et al. 2003). 



2 UNIVERSAL DARK MATTER DENSITY 
PROFILE 

The dark matter density profile, Pdm{f) suggested by many 
high resolution N - body simulations is well described by a 
self-similar form. We assume that the gas mass is negligible 
compared to the total dark matter mass and adopt this uni- 
versal density profile for dark matter in clusters. The profile 
is expressed in terms of a characteristic radius (e.g., in 
Komatsu & Seljak 2002): 



Pdm = psydm{r/rs) 



(1) 



where pa is a normalization factor which represents a char- 
acteristic density at a characteristic radius, r — rs- This 
characteristic radius describes a typical scale at which the 
profile slope changes from the outer value to the inner value. 
The functional form of yam (x) is given by 



Vdmix) 



(2) 



a;«(l -|-a;)3-« 

Here the parameter a characterizes the shape of the 
profile. Since the dark matter density profile is self-similar, 



the dark matter mass profile is also self-similar. So, the dark 
matter mass enclosed within a radius r is 

M(< r) = 47rpsr?m(r/rs) (3) 

where, m{x) is a non-dimensional mass profile given by 



rx 
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m{x) = / dui^ydm{u) = ln(l -|- x) 



{i + xY 



(4) 



Here, the last equality is valid for a — 1 which is the 
much used NFW profile (Navarro et al. 1996, 1997), the in- 
tegral being evaluated by Suto et al. (1998). 

The definition of the virial radius, Rvir, is the radius 



Myir = M{<c), where 



(5) 



is a dimensionless parameter called the 'concentration pa- 
rameter'. Evaluating equation (3) at the virial radius, the 
normalization factor, ps, is fixed at; 



Ps= C 



(6) 



The virial radius, Rvir {M„ir,z) is calculated with the 
spherical collapse model (Peebles 1980), 



Rvi' 
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(7) 



where the second equality comes from evaluating Rvir from 
equation (6). Here Ac{z) is the spherical over-density of the 
virialized halo within Rvir at in units of the critical den- 
sity of the universe at z, Pc{z). Following Komatsu & Seljak 
(2002), we assume a value Ac(z = 0) = 100 for a cosmolog- 
ical model with Q.m ~ 0.29 and JIa = 0.71. 

We follow Bullock et al. (2001) in adopting the approx- 
imation for c as a function of the virial mass of the cluster. 
They give the median values of 'c' and also the la devia- 
tions: 



k(- 
Vi. 



5 X 10i3/i-iMp 



(8) 



with K = 9 reproducing the best-fit and K = 13.5 and 
K = 5.8 reproducing the -|-lcr and the —la values in the 
concentration parameter. These values of the concentration 
parameter arc also consistent with the findings of Seljak & 
Huffenberger (2003). 

The above set of equations specify the dark matter den- 
sity profile of a particular mass cluster. Next, we turn our 
attention to the density profile of the gas in hydrostatic equi- 
librium with this dark matter distribution. 

To compare our results with observations, which usu- 
ally uses the radius R200 where the over-density is 200, we 
compute this radius in each case and present our results in 
the terms of 7i2oo- 



3 HYDROSTATIC EQUILIBRIUM OF GAS 
AND DARK MATTER 

Our aim in this section is to determine the density profile of 
the intra-cluster gas using the universal temperature profile 
(discussed later) and assuming that the gas is in hydrostatic 
equilibrium with the background dark matter potential. 
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The typically smooth morphology of the X-ray emission 
from the hot intra-cluster medium leads naturally to the 
hypothesis that the gas is near equilibrium, stratified along 
isopotential surfaces in a mildly evolving distribution of dark 
matter, gas and galaxies. This suggests that the assumption 
of hydrostatic equilibrium for such relaxed clusters is mostly 
justified. 

3.1 Universal Temperature Profile of Gas 

The "universal temperature profile" used for our calculation 
(Loken et al. 2002) is (normalized by the emission-weighted 
temperature): 

T To 



(T> (1 + 



(9) 



where (T) is the emission-weighted temperature of the clus- 
ter. To = 1.33, Ox = Rvir/1-5, and 5 = 1.6 on the radial 
range (0.04-1.0) Ryir- To determine the emission- weighted 
temperature from the cluster mass, we use a relation that 
arises from adiabatic evolution of the gas in cluster. Af- 
shordi & Con (2002) liave shown that the observations of 
Finoguenov ot al. (2001) of M5oo~{T) relation in clusters 
can be understood from gravitational processes alone. We 
therefore use this empirical relation (M5oo-(T')) derived by 
Finoguenov et al. 2001: 



(10) 



M5oo = (2.64l°-)10-Me(^)'"-°°^ 

where kb is the Boltzmann constant and M500 has been cal- 
culated self-consistcntly by taking the total mass within the 
radius where the over-density is 5 > 500. 

The main motivation of using this universal tempera- 
ture profile is to see if this profile which arises just out of 
gravitational interactions alone in clusters predicts anything 
different from the previously used default temperature pro- 
files from gravitational interactions alone. This temperature 
profile is a result of a high resolution simulation, without any 
input from non-gravitational processes, carried out by Lo- 
ken et al. (2002) which makes use of a Eulerian-beised, adap- 
tive mesli-rofinomont code that captures the shocks that are 
essential for correctly modelling cluster temperatures (for 
details of the simulation, refer to Loken et al. (2002)). The 
temperature profiles of the simulated ACDM and SCDM 
clusters are remarkably similar and are well fit by this uni- 
versal temperature profile. 

Finoguenov et al. 2001 also point out that the M5oo-(T) 
relation becomes flatter (M500 oc T^'^®) when low mass clus- 
ters (M500 < 5 X 10" Mo) arc excluded i.e. the Mboo-T re- 
lation becomes closer to self-similar relation of M500 oc T^''\ 
We have tried out with a flatter Msoo-T relation and seen 
that the results do not change appreciably. The code used 
by Loken et al. (2002) was tested as a part of the Santa Bar- 
bara cluster comparison project (Fronk et al. 1999) in which 
12 groups simulated a Coma-like cluster using a variety of 
codes and resolutions. They have shown that their results 
are among the highest-resolution results presented in the par- 
per(central resolution of 7.8h~^kpc). They have also pointed 
out that their results for the Santa Barbara cluster arc in 
excellent agreement with the results obtained from those of 
a new, completely independent code (Kravtsov, Klypin & 
Hoffman 2002). 



This profile is in good agreement with the observational 
results of Markevitch et al. 1998 but diverges, primarily in 
the innermost regions, from their fit which assumes a poly- 
tropic equation of state. This profile is also in very good 
agreement with a recent sample of 21 clusters observed by 
BeppoSAX (De Grandi & Molendi 2002) with and without 
cooling flows. Although the simulation result is consistent 
with the data at outer radii, there is some difference in the 
inner region of clusters, indicating that there could be some 
additional physics at small radii (r < O.l-R^ir) (Nath 2003). 

We emphasize that this temperature profile does not 
include the effects of cooling and galaxy feedback or for that 
matter any additional physics. This temperature profile can 
be used, therefore, to probe the influence of gravitational 
interactions on the ICM, at radii > 0.17?«ir. 

It is instructive to compare the above mentioned tem- 
perature proflle with the temperature profile calculated by 
assuming that pgasir) = fgasPdm(r), Vr where fgas = 0.105 
(e.g in Bryan 2000, Wu & Xue 2002b) for a range of cluster 
masses. This is the self-similar model that has been used 
as a calibrator for the influence of gravitational processes. 
Although one does not expect in reality for the above pro- 
portionality to hold for arbitrarily small radii, it has boon 
expected that the influences of shocks resulting from gravita- 
tional interactions alone does not extend beyond ~ O.lRvir 
(e.g., Bryan 2000). This model therefore has been widely 
used to calculate the expected entropy of the gas at radii 
> O.lRvir from gravitational interactions alone, and com- 
pare this expectations with the observed values. 

We compare the temperature profiles assumed for the 
present model, and the profiles obtained using the self- 
similar model {pgas{r) oc pdm{r)) in Figure 1, normalizing 
the temperature profiles by T200, where 



7200 



GM-zaaj-irrip 



2R20 



(11) 



where, M200 = /g 4-Kpdmr^dr. 

It is scon that the universal temperature profiles (see 
Figure 1) flatten towards the inner regions of the cluster 
(within 0.2i?20o), whereas the temperature proflles calcu- 
lated using pgas(r) oc pdrn(r) (hereafter, referred to as self- 
similar profiles) dip at the central region of the cluster. 



3.2 Density Proflle of Gas 

We next numerically evaluate the density profile of the intrar 
cluster gas with the above defined temperature profile given 
in equation (9) assuming hydrostatic equilibrium. 

Let us consider a spherical gas cloud with temperature 
T(r); then its density distribution pg{r) in hydrostatic equi- 
librium satisfies 



d 



GM{< r) 



(12) 



(13) 



pgas{r) dr 
where, 

'pgasir)' 

fj,mp 

where M(< r) is the total mass inside radius r (equation 
3) and p and mp denote the mean molecular weight (we 
use fj, = 0.59) and the proton mass. The boundary condi- 
tion which is imposed in evaluating this integral is that the 



Pgas = (^^^)kbT{r) 
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Figure 1. The variation of gas temperature (scaled by T200) with 
scaled radius for 2 clusters of different masses. The region shaded 
with vertical lines is the spread in gas temperature profiles of Lo- 
ken ct al. (2002) due to la spread in the concentration parameter 
'c' for a cluster of (T) 8.5 keV with the solid line in the middle 
being the profile for the best fit value of 'c'. The region shaded 
with solid slanted lines is the universal temperature profile with 
the same spread in 'c' for a poor cluster of (T) 0.85 keV with 
the long-dashed line in the middle being for the best fit value of 
'c'. The region shaded dark with closely spaced horizontal lines is 
the result of the self-similar calculations for a cluster (T) 8.5 keV 
representing the spread in gas temperature for a la spread in 'c' 
with the solid line in the middle being the result for the best fit 
'c'. Finally the region shaded with broken dashed lines is the gas 
temperature for tlie self-similar model of a cluster of (T) 0.85 keV 
with the spread being due to a spread in 'c'. 

gas-fraction, fgas , within the virial radius is universal and is 
equal to 0.105, as recently found by Ettori (2003) for a sam- 
ple of low and high redshift clusters. Since the total gas mass 
is negligible compared to the dark matter, Mtotai ^ Mdm, 
mass in dark matter, and therefore fgas ~ T^f^ ■ "^^^ ques- 
tion whether or not fgas is independent of cluster mass has 
been a topic of debate in the literature (White & Fabian, 
1995; Ettori & Fabian, 1999; Markevitch et al. 1999; Wu & 
Xue, 2000). Although recent observations by Sanderson et 
al. (2003) show an apparent trend of fgas being smaller for 
lower temperature systems, they also found that a universal 
value of fgas can fit their data. According to them, the un- 
weighted mean of their data set of gas-fraction within the 
virial radius is a constant close to fgas = 0.13 ± 0.01 to 
0.1 ± 0.01, the variation being due to the variation in the 
slope of the M-{T) relation. 

We note that this normalization provides a conservative 
estimate of the entropy of the gas from gravitational collapse 
alone since a lower value would only increase the entropy at 
all radii. 

Figure 2 shows the gas-density profiles of 2 clusters 
of different masses. It compares the density profiles of the 
present model (lower set of curves) to the density pro- 




0.01 0.1 1 

Scaled Radius (R/Ra^,) 

Figure 2. Gas density profiles for 2 different clusters with scaled 
radius. The region shaded with horizontal lines is the spread in gas 
density as a result of the la spread in the concentration parameter 
for a rich cluster of (T) = 8.5 keV. The region shaded with slanted 
lines corresponds to a low mass cluster (0.85 keV). The higher set 
of curves correspond to the self-similar model (pgas = fgasPdm)- 
The dot-dot-dashed line represents the cooler cluster (0.85 keV) 
and the dashed line the hotter of the two (8.5 keV). 



files from self-similar assumptions (higher set of curves) 
with a constant proportionality factor of fgas = 0.105 i.e., 
pgas = fgasPdm- As expcctcd, it is sccu that the gas density 
is much shallower in the inner parts of the cluster as a result 
of the universal temperature profile being flat at the inner 
regions. Interestingly enough, the density profiles obtained 
from the universal temperature profile deviates from self- 
similar expectations at radii much larger than 0.1-R2()()- For 
a poor cluster with emission weighted temperature of 0.85 
keV the deviations become significant even at r ~ 0.4i?2oo ■ 

We note that the emergence of a core in the gas density 
distribution has been noticed by previous authors of numer- 
ical simulations, even in the absence of non-gravitational 
heating and cooling processes (Frenk et al. 1999). This ap- 
pears to result from the transfer of energy between baryonic 
and dark matter during merger events (Eke et al. 1998). 

We have also done all the calculations described here 
for a coma- like cluster {M^ir = 1.1 x W^^Mq) for a S-CDM 
universe using the same parameters as in the Santa Barbara 
Cluster Comparison Project (Frenk et al. 1999) to check for 
the consistency of our method and results. We find that all 
the properties like gas-density profile. X-ray luminosity cal- 
culated using the prescription described in this paper match 
the simulated results of Frenk et al. (1999). 
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4 IMPLICATIONS OF THE UNIVERSAL 
TEMPERATURE PROFILE AND THE 
DERIVED GAS DENSITY PROFILE 

In this section, we focus on the imphcations of the above 
temperature and gas-density profiles on the other physical 
properties of this intra-clustor gas like entropy, gas-mass, 
variation of Mgas with cluster mass or emission-weighted 
temperature, T and gas-fraction, fgas in light of recent ob- 
servations. 

4.1 Entropy Profiles and Scaling Properties 

For convenience, 'entropy' for the intra-cluster gas is defined 
as 

S=^(keVcm^) (14) 

Tie 

where, T is the temperature of the gas and rie the parti- 
cle density. This quantity is directly related to observations. 
This has been referred to by a number of authors as 'adia- 

bat', since (apart from a constant relating to mean particle 
mass) it is the coefficient relating pressure and density in 
the adiabatic relationship P = Kp^ . Hence S is conserved 
in any adiabatic process. Note that the true thermodynamic 
entropy is related to this definition via a logarithm and ad- 
ditive constant. 

In this section, we discuss the scaled entropy profiles 
(scaled with the emission- weighted temperature, (T)) ob- 
tained from our density and temperature profiles for 5 dif- 
ferent mass clusters. We calculate the emission-weighted 
temperature corresponding to the profiles discussed above 
within a fiducial radius of 0.3i?2oo (in the band 0.5 — 10 
keV), using the Raymond Smith code, for a metallicity of 
Z/Zq = 0.3: 

^ J°"''''Unr'ni{r)n.{r)eo.,-roT{r)dr 

A7rr^ni(r)ne(r)eo.5-iodr 

where rii and Ue represent the ion and electron density and 
eo.5-10 denotes the emissivity relevant for the 0.5 — 10 keV 
band. We find that the emission weighted temperature ob- 
tained in this manner matches well (within 0.5%) with the 
value assumed to calculate the temperature profile itself 
(from equation 9). This shows that the system of equations 
used for our calculations is self-consistent. 

Under the assumption that all these systems form at the 
same redshift, their mean mass densities should be identical. 
Hence in the simple self-similar case, where all have similar 
profiles and identical gas-fractions, S will simply scale with 
emission- weighted temperature (T). We apply this scaling 
and scale the radial coordinate to R200 for each system, de- 
rived as mentioned above. 

We show these scaled entropy profiles in Figure 3. It 
can be seen that the entropy profiles of the cooler systems, 
scaled in the above mentioned way tend to be significantly 
and systematically higher than that of rich clusters. In these 
derived entropy profiles, we notice that they generally flatten 
in the very interior parts of the clusters (inside 0.05J?2oo) 
resulting from the flattening of density distribution at these 
radii. It is also seen that there is a noticeable tendency for 
the scaled entropy to be higher, at a given scaled radius, in 
cooler systems. Simulations and analytical models of cluster 



1000 r 1 1 1 1 1 — I — I — r-| i 1 1 1 1 — I — I — r 




Scaled Radius (R/Rj,,,,) 



Figure 3. Scaled entropy profiles (scaled by I/T200) for 2 clus- 
ters. The region shaded with horizontal lines represents the spread 
in the scaled entropy profiles due to the Isigma spread in the 
concentration parameter 'c' for a cluster with (T) of 8.5 kcV and 
the region shaded with slanted lines represents the spread in the 
scaled entropy profiles for a cluster with (T) of 2.44 keV. The 
region shaded with vertical lines is for a low mass cluster of (T) 
0.85 keV. The bottom line (solid) shows the slope of 1.1 expected 
from shock heating. Its normalization is arbitrary. 



formation involving heating from accretion shocks, produce 
entropy profiles with logarithmic slopes of approximately 
1.1 (Tozzi & Norman 2001), which agrees rather well with 
the slope of the calculated profiles outside R « 0.2_R2oo, for 
rich clusters but it does not show good agreement with the 
entropy profiles for poor clusters. 

The general trend of our calculated scaled entropy pro- 
files are in good agreement with the recent results of PSF03. 
However, it is seen that the values in general are systemati- 
cally lower than that of PSF03 at around 0.01J?2oo and also 
at the outer reaches of the cluster. 

In Figure 4, we plot the variation of entropy S at 
0.1i?2()() with emission- weighted temperature, (T). The data 
points with the error bars are from PSF03. The region 
shaded with oblique lines represents the entropy calculated 
from the above described model with the la spread in the 
concentration parameter "c' and the region shaded with 
crossed lines corresponds to the self-similar density profiles. 

The discovery of an entropy fioor in galaxy groups and 
clusters (Ponman et al. 1999) was based on the measure- 
ment of gas entropy at 0.17?2oo, in systems spanning a wide 
temperature range. This radius was chosen to lie close to 
the centre, where accretion shock-generated entropy should 
be minimum, hence maximising the sensitivity to any ad- 
ditional entropy, whilst lying outside the region where the 
cooling time is less than the age of the universe, and hence 
the entropy may be reduced. This initial study was improved 
by Lloyd-Davics ct al. (2000) who derived an entropy floor 
value of ISdhgg^^keVcm'^ from a sample of 20 systems. How- 
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Temperature (keV) 

Figure 4. Gas entropy at 0.1i?2()() as a function of emission- 
weighted temperature. The region shaded with slanted lines shows 
the spread in entropy due to the la spread in the concentration 
parameter with the solid line in the centre being the results of 
the best-fit 'c' for the above described model, the data points 
are from PSF03 and the region shaded with crossed lines shows 
the spread in entropy due to the spread in 'c' for the self-similaj: 
model {pgas oc Pdrn ) with the dashed line in the centre being the 
results of the best-fit 'c'. 
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Temperature (keV) 

Figure 5. Entropy at R-,oo as a function of emission- weighted 
temperature. The solid line is the result of the above described 
model with the best-fit value of 'c' and the region shaded with 
slanted lines being the spread in the entropy due to the la spread 
in 'c'. The data points are from PSF03. The dashed line is the 
entropy calculated from the self-similar model(/Ogas oc Pdm) with 
the best-fit value of 'c' and the region shaded with crossed lines 
showing the spread in entropy due to the spread in 'c'. 



ever, the recent results of PSF03 show that there is no such 
entropy floor. They point out that an unweighted orthog- 
onal fit to the data points, which have been grouped into 
temperature bins, gives a logarithmic slope of 0.57 ± 0.04, 
as opposed to 5 oc T in self-similar relation. 

As is seen in the figure, the entropy calculated from our 
model is higher than the previously calculated entropy from 
self-similar models (Wu & Xue, 2002b) for the poor clusters 
and similar or slightly lower for rich clusters. However it is 
lower than the data points from PSF03, the difference being 
more pronounced for poor clusters. This is because of the 
fact that the density profile deviates from the self-similar 
models even at O.lRvir. 

We have noted earlier that the gas density profiles cor- 
responding to the universal temperature profile is substan- 
tially flatter than the self-similar models even at radii larger 
than 0.1ii2oo (see Figure 2). It is therefore instructive to 
compare the entropy at larger radii with that from self- 
similar models and data. Recently Finoguenov et al. (2002) 
compared the entropy at R500 (which is ~ |il200 for the 
range of cluster masses) with those expected from self- 
similar models and concluded that there is excess entropy 
even at this large radius, indicating the large scale influence 
of preheating processes. 

In Figure 5, we show a plot of 'entropy' S{R5oo) with 
emission- weighted temperature, (T). The region shaded 
with slanted lines represents the present model with the la 
spread in the concentration parameter and the region shaded 
with crossed lines shows the entropy calculated from the self- 
similar model {pgas{r) oc pdmir))- The solid line through the 
middle of the region shaded with slanted lines represents the 
results for the best- fit 'c'. Data from PSF03 arc also plot- 
ted for comparison. We find that the level of entropy at 
R500 is reasonably consistent with the observed values for 
rich clusters but they are lower than the observed values for 



intermediate and low mass clusters, with the deficiency bo- 
coming appreciable for poor clusters. We note that previous 
authors had re-normalized the expectations for self-similar 
models by matching them with the entropy of the richest 
clusters and thus concluded the presence of excess entropy 
at i?5()(). Wc do not normalize our calculated entropies in 
this manner in this paper. However, it still shows that there 
is a need for non-gravitational heating even at large radii 
especially for low mas clusters to fulfill the requirement of 
this excess entropy. 

The Figures 3 and 4 clearly show that the entropy of 
the ICM in the central regions (< 0.1i?2oo) from gravita- 
tional processes alone, is larger than the previous expec- 
tations from self-similar models. It is then reasonable to 
conclude that the entropy imparted by interactions between 
dark matter and baryons (Eke et al. 1998) has been under- 
estimated. Therefore, we find that gravitational interactions 
impart more entropy to the gas in the central regions than 
estimated earlier but the problem of excess entropy for low 
mass clusters still remains. It is interesting to find that the 
entropy at JJsoo in the present model is consistent with the 
observed values for rich clusters but are low for poorer clus- 
ters which probably confirms the requirement of some pro- 
cess which would help preferentially increase their entropy. 
It is not surprising to find consistency in the case of rich clus- 
ters though, as we have already noted that the temperature 
profile of Loken et al. (2002) is consistent with the observed 
profiles at outer radii. 



4.2 Gas Distribution 

In this section, we discuss the effects of this temperature 
profile and the resulting density profile on the relation be- 
tween Mgas with mean emission-weighted temperature and 
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Temperature (keV) 

Figure 7. Gas fraction within O.3R200 a-s a function of emission- weighted temperature, (T). The region shaded with slanted Unes 
represent the above discussed model with spread in 'c' with the solid line in the centre being the results of the best-fit value of 'c' and 
the dashed line represents the model with pgas oc Pdm- The data points are from Sanderson et al. 2003. 




T (keV) 

Figure 6. Comparison of Mgas (-Rsoo) — T relations. The data 
points with error bars represent gas mass determinations of 
MME99 within R500. The solid line is the result of the present 
model using the best-fit value of 'c' with the region shaded with 
slanted lines representing the spread in 'c', the region shaded with 
crossed lines being the prediction of the self-similar model with 
the dashed line in the centre being the results of the best-fit value 
of 'c'. 



the gas-fraction, fgas profiles and the variation of fgas, in 
the inner regions, with emission-weighted temperatures. 



4.2.1 Mgas (Rboo) - (T) 

In figure 6, we present the Mgas — (T) relation as predicted 
by the present model (solid line), the relation derived from 



the self-similar model (dotted line) , and the data points from 
Mohr et al. (1999) within -Rsoo. It is seen that the gas mass 
within i?5oo calculated from the present model is slightly 
higher than the data points, but lower than the expectations 
from self-similar models for clusters with (T) < 3 keV and 
slightly higher for clusters with (T) > 3 keV, which was 
previously thought to be the result of gravitational processes 
alone. The logarithmic slope of our curve (~ 2.08)is steeper 
than the self-similar slope of 1.5, and close to the observed 
slope. 



4.2.2 Gas fraction fgas and its spatial variation 

The variation in the gas fraction (fgas), evaluated within a 
characteristic radius of 0.3i?2oo is shown is Figure 7. The 
solid line is the calculated gas fraction from the present 
model and the dotted line shows the results of self-similar 
model. The data have been taken from PSF03. It is noted 
here that the gas fraction obtained from the present model 
is slightly higher than the data points. However, it is also 
seen here that there is a clear trend for cooler systems to 
have a smaller mass fraction of X-ray emitting gas in the 
central regions, as a result of gravitational processes alone. 
The gas fraction obtained from the self- similar model is a 
constant as it should be by definition. 

To understand the behaviour of gas fraction with ra- 
dius better, we have plotted fgas{r) with the scaled radius, 
r/R'ioo in Figure 8 for 5 clusters of difiierent masses. The 
general trend seen here is for gas-fraction to rise monotoni- 
cally with radius (especially for clusters with (T) < 3 keV) 
all the way till -R2oo. There is, however, a flattening of the 
profiles for the rich clusters ((T) > 3 keV) beyond 0.5i?2oo 
with a slight bump around 0.3J?2oo. It should be noted that 
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Figure 8. Spatial variation of gas fraction witliin a given radius 
(normalized to i?20o) for three different clusters with different 
emission-weighted temperature T. The solid line with the region 
shaded with horizontal lines represents the coolest system (0.85 
keV), increasing in temperature through long-dashed with a re- 
gion shaded with slanted lines (2.22 keV) and finally short-dashed 
(8.5 keV) line with the region shaded with vertical lines. These 
shaded regions represent the spread in gas fraction due to the Icr 
spread in 'c'. 



the fit provided by Lokcn ot al. 2002 is accurate to about 10 
% for r < 0.5Rvir and underestimates tlieir simulated tem- 
perature profile in this region. This may increase the gas 
density profile at smaller radii and account for the bump in 
Figure 8. It can be seen clearly that the profiles lie in or- 
der of temperature such that, at a fixed radius, gas fraction 
decreases as temperature decreases, mirroring the trend in 
Figure 7. 

If cluster evolution (with gravitational processes alone) 
were an entirely self-similar process, then fgas should be a 
constant at a given over- density in all objects. Figure 9 
clearly shows that this is not the case, even for evolution 
in the presence of gravity alone. Poor clusters seem to have 
a much lower gas mass fraction compared to rich clusters 
at the same over-density which was also seen by David et 
al. 1995. Again, a bump is seen for the richest cluster as seen 
in Figure 8. We have already discussed the probable reasons 
for this feature earlier. 



4.3 X-ray Luminosity- Temperature Relation 

In this section, we compute the bolometric X-ray luminosity, 
corresponding to the profiles discussed above (in the band 
0.5-10 keV), using the Raymond-Smith code, for a metallic- 
ity of Z/Zq = 0.3. We compute the luminosities within the 
virial radius 0.3J?2oo- The X-ray luminosity is computed as, 

/•0.3-R200 

Lx = 47rr^ni(r)ne(r)eo.5-io (16) 

Jo 




10' 1000 
s 



Figure 9. Gas mass fraction as a function of the over-density S, 
for five different clusters. The solid line represents the coolest sys- 
tem (0.81 keV), increasing in temperature through dashed (1.26 
keV), dot-dashed (2.22 keV), dotted (3.29 keV) and finally dash- 
dot-dot-dotted (7.36 keV) line. 



where rii, He represent the ion and electron density and 
eo.5-10 denotes the emissivity relevant for 0.5-10 keV band. 
We present the results in Figure 10. 

It is seen, from Figure 10, that the luminosity calcu- 
lated from the present model (with the above quoted values 
of the cosmological parameters) is close to the observed data 
(scaled with h = 0.71 as used for the models) for clusters 
with emission- weighted temperature (T) above 1.0 keV. The 
data points are from Arnaud & Evrard 1999 and Markevitch 
1998 for clusters above 5.0 keV and from Heldson & Pon- 
man 2000 for low (T) regime i.e. for low mass clusters and 
groups are also shown. Interestingly, compared to the dotted 
line which represents the self-similar model {pgas oc pdm), 
the present model shows that the luminosity for low mass 
clusters (below 5 keV) is closer to data points. However, the 
luminosity is still somewhat over-estimated in this model 
when compared to data. This again indicates that there is 
requirement for some non-gravitational heating (preferen- 
tially in low mass clusters) to reduce the gas density further 
and thus the X-ray luminosity. We note here that tiie X-ray 
luminosity depends strongly on the assumed metallicity, es- 
pecially for low temperature systems, and the uncertainty 
over abundance of gas in poor clusters is yet to be resolved 
(Buote 2000, Davis 1999). 

In Figure 11, we have plotted the X-ray luminosity 
integrated within a radius iiiooo where the over-density 
6 > 1000. The results of our calculations are compared with 
the best-fit results of Ettori et al. 2002. It is interesting to 
find that the X-ray luminosities calculated from the above 
model lie within the la spread of the results of Ettori et 
al. 2002. 
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T (keV) 

Figure 10. Relation between bolometric X-ray luminosity Lx 
and emission- weighted temperature ({T}). The data points rep- 
resented by 'stars' show measurements of clusters with insignif- 
icant cooling flows compiled by Arnaud & Evrard (1999). Open 
squares show cooling flow-corrected measurements by Markevitch 
et al. (1998). The data points with error bars show group data 
from Helsdon & Ponman (2000). The shaded region represents 
X-ray luminosity calculated using the above model with the la 
spread in 'c' with the solid line representing the median value and 
the region shaded with crossed lines is the result of the self-similar 
model with the same spread in 'c'. The models assume a ACDM 
cosmology with Qm = 0.29, Qa = 0.71, and Q^, = 0.047, and a 
Hubble parameter of /i = 0.71 has been applied to the models 
and the data. 




T (keV) 

Figure 11. Relation between bolometric X-ray luminosity Lx 
and emission- weighted temperature {(T)). The shaded region 
enclosed in between the two solid lines represents the Icr spread 
which are the results of the best-fit analysis of the data of Ettori 
et al. 2002 with the solid line through the middle being the results 
of using the best fit parameters. The region shaded with slanted 
lines is the la spread calculated using the la deviations in the 
concentration parameter stated above with the dashed line being 
the results of the median value of 'c' and the region shaded with 
crossed lines in the result of similar calculations using the self- 
similar model. 



5 DISCUSSION 

To recapitulate, our attempt here has been to study in de- 
tail the implications of the universal temperature profile ob- 
tained from recent high resolution simulations (Loken et 
al. 2002), with no input from non-gravitational processes, 
such as heating or cooling. Wc have compared the predic- 
tions of entropy and gas fraction from this temperature pro- 
file with the self-similar model {pgas oc Pdm), which has been 
used in the literature as a reference model for calibrating the 
infiuence of gravitational processes, at radii r > O.lRvir- The 
only assumption that wc have made in this paper is that of 
hydrostatic equilibrium, apart from the assumption of the 
background dark matter density profile. We have found that 
the result of gravitational processes alone is much different 
from what has been used in the literature so far. 

Firstly, the widely used assumption that entropy im- 
parted to the gas from the interaction between dark mat- 
ter and baryons is limited to the very central region, r <C 
O.lRvir appears to be simplistic. The corresponding assump- 
tion that the gas density profile expected from gravitational 
processes alone is proportional to the dark matter density 
profile at r > O.lRyir seems to be violated. Curves in Figure 
2 show that the gas density profile deviates from the simple 
proportionality at a much larger radii, r > 0.5i?2oo, and this 
deviation is larger for lower temperature systems. 

If the temperature profile obtained by Loken et 
al. (2002) is confirmed to be the one that is expected from 
gravitational processes alone, then the implications dis- 
cussed here are inevitable. Instead of using the self-similar 
models, one must use the present model to benchmark the 
expectations from gravitational processes, and then compare 
it with the data to infer the need for any additional physics. 

However, to conclude definitively about the presence of 
any non-gravitational heat in the ICM, one needs to know 
better about the dark matter profile in clusters and to re- 
solve the uncertainties in the concentration parameter in 
dark matter halos of clusters. This idea was also explored 
by Lloyd-Davies et al. (2002) who concluded that the scatter 
in the concentration parameter plays a role in the determi- 
nation of the contribution of non-gravitational heating in the 
intra-cluster medium. Also, in this context, it is worthwhile 
to point out that the dependence of 'c' on redshift would 
have an efi^ect on the results. This is because the assump- 
tion that the clusters were all formed at a redshift «/ = is a 
simplification of the actual picture. Different clusters would 
form at difi'erent redsliifts according to the theory of struc- 
ture formation (Press & Schechter, 1974) and that would 
mean that the concentration parameters would be different. 
The values of 'c' would be lower as c oc (as found 

by Bullock ct al. 2001). This will make the results slightly 
different, e.g., entropy at 0.1ii2oo would be higher than the 
values plotted now and thus closer to the data and the X- 
ray luminosity would also be lower and thus closer to the 
data points. liowcvor, the entropy at 7?5oo will increase as a 
result of this and the gas fraction at 0.37?200 will decrease, 
thus being closer to the respective data points. 

We find that for entropy at 0.1J?2oo the present model 
produces lower entropy than the data for all clusters. The 
difference is larger for poorer clusters. The problem would 
then be to increase the entropy, especially for poor clusters, 
as has been required from theoretical models earlier. The 
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requirement of non-gravitational heating as estimated from 
earlier theoretical models would be reduced because as it 
is seen here that the entropy at inner regions were under- 
estimated in earlier models. We also find that the entropy 
expected at RrMo is consistent with the data for rich clusters 
but it is lower than the data points for low mass clusters em- 
phasizing the need for some process to heat the gas even at 
large radii. It is however interesting that the scaled entropy 
profiles are similar to that observed (PSF03, Mushotzky et 
al. 2003). 

It is possible that the discrepancy between the present 
model and the data for entropy at 0.1J?2oo and X-ray lumi- 
nosity for rich clusters can be alleviated by gaseous processes 
such as thermal conduction. As Lokcn ct al. (2002) com- 
mented, the observed temperature profile is more flattened 
than their simulated profile, with somewhat larger tempera- 
ture at ~ 0.1i?2oo. This could be due to thermal conduction 
(e.g., Nath 2003). This would then decrease the density in 
the inner regions and would (1) increase the entropy & (2) 
decrease X-ray luminosity to be consistent with data for 
rich clusters. We note that thermal conduction is however 
less important poor clusters. 

We have noted that normalizing the gas content of clus- 
ters with a constant gas fraction for all clusters provides a 
conservative estimate of entropy from gravitational collapse 
alone. Also, we have found that changing the exponent of 
the M500-T relation to 1.5 does not change the results much. 

The total gas mass at a fiducial radius of -Rsoo expected 
from the present model are again closer to the data than the 
previous self-similar model. The gas fraction at 0.3i?2oo cal- 
culated from this model is higher than the data but it agrees 
better in comparison to the previously calculated ones from 
the self-similar model. It also rises with the temperature of 
the cluster as observed. 

There is, therefore, some difference between the ex- 
pectations from the present model and the data. If the 
present model is a realistic indicator of gravitational pro- 
cesses alone, then the results from this model should be 
compared with the data, to determine the requirements of 
additional physics, if any, to explain the data. 



6 CONCLUSION 

The primary aim of this work was to study the implications 
of the 'universal temperature profile' arising out of pure 
gravitational interactions in gala^xy clusters. We have deter- 
mined the gas density profile corresponding to this tempera- 
ture profile, and studied various implications of this profile. 
The only assumptions made in this paper is that of hydro- 
static equilibrium and the temperature profile of Loken et 
al. 2002. 

We have also shown the dependence of the above results 
on the uncertainty in the knowledge of the concentration 

parameter 'c'. 

Given the uncertainty in the concentration parameter, 
the main results are summarized below: 

(a) Gas density profiles expected from gravitational pro- 
cesses alone is flatter than previously thought, even at radii 
much larger than O.lRyir. 

(b) Entropy expected firom gravitational processes alone 
at 0.1ii2oo are larger than previously thought, especially for 



low mass clusters, but still lower than the observed values, 
with the discrepancy increasing for low mass clusters. The 
entropy expected at R500 is consistent with observed values 
for rich clusters but it is lower than the data points for low 
mass clusters. Thus it emphasizes the need for non- gravi- 
tational heating preferentially for low mass clusters even at 
large radii. 

(c) Gas fraction in the inner region (0.3i?2oo) expected 
from gravitational process alone is much smaller than pre- 
viously thought, and but slightly higher than the observed 
values. 

We therefore conclude that if the temperature profile 

of Loken et al. (2002) is indeed the result of evolution of 
the intracluster gas involving gravitational processes alone, 
then the contribution of non-gravitational processes to the 
physics of ICM would have to be revised. 
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